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A new method to drive a system of neutral dipolar fermions into the lowest Landau level regime 
is proposed. By employing adiabatic spin-flip processes in combination with a diabatic transfer, the 
fermions are pumped to higher orbital angular momentum states in a repeated scheme that allows 
for the precise control over the final angular momentum. A simple analytical model is derived to 
quantify the transfer and compare the approach to rapidly rotating systems. Numerical simulations 
of the transfer process have been performed for small, interacting systems. 
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Despite being ideal models for complicated solid state 
systems, ultracold quantum gases lack one important as- 
pect of the electronic complex: because of the charge neu- 
trality of the atoms, there are no mobile charge carriers 
that possess a direct coupling to the magnetic vector po- 
tential. Plenty of interesting effects, however, arise when 
charged particles are subject to high magnetic fields in 
low dimensional systems. The most prominent ones are 
the integer quantum Hall effect [1] as an example for the 
appearance of topological states, as well as the fractional 
quantum Hall effect [2], potentially giving rise to funda- 
mental excitations with non-Abelian statistics. 

Several schemes have been proposed to simulate the 
effect of magnetic fields for neutral particles. Artificial 
gauge fields can be created by imprinting phases, mak- 
ing use of the Peierls substitution in optical lattices [3- 
5], or by tailoring spatially dependent Hamiltonians to 
generate geometric phases [6]; for an overview see [7]. 
Rapidly rotating quantum gases provide an alternative 
route via Larmor's theorem [8, 9], Several theoretical 
proposals demonstrate the appearance of highly corre- 
lated quantum Hall states for dipolar bosons [10] and 
fermions [11, 12]. However, the experimental realization 
of quantum Hall states has been elusive so far. For ro- 
tating systems, the main problem is the precise control 
on the rotation frequency, which is required to reach the 
lowest Landau level without crossing the rotational in- 
stability [13]. 

Here we propose a new scheme to access the regime 
of fast rotation for a dipolar Fermi gas such as 161 Dy, 
which has recently been cooled to the quantum degen- 
erate regime [14]. Starting from a spin-polarized state, 
dipolar interactions can lead to spin relaxation with a net 
angular momentum transfer [15]. This is known as the 
Einstein-de Haas effect [16] and has been proposed to 
create rotating Bose-Einstein condensates [17, 18]. We 
suggest using this mechanism in a trapped, quasi-two- 
dimensional system to control the amount of angular 
momentum, and - by repeated application of the trans- 
fer scheme - reach the lowest Landau level (LLL). This 



scheme allows for direct control over the total angular 
momentum instead of the rotation frequency and cir- 
cumvents the prime experimental difficulties toward the 
realization of the quantum Hall regime in harmonically 
trapped gases. 

We consider a system of N fermionic atoms with mag- 
netic dipole moments fi. While extensions to schemes 
with polar molecules are possible, the permanent dipole 
moments of the atoms lead to some simplifications. To 
shorten the discussion, we consider only two internal lev- 
els (pseudospin 1/2). The particles are confined in a 
quasi-two-dimensional harmonic trap with a radial fre- 
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FIG. 1. (a) Dipolar particles, trapped in a quasi-2D geom- 
etry with a radial confinement ui. When the external mag- 
netic field B is tuned in resonance, dipolar interactions Vdd 
can induce spin relaxation processes, leading to a net angular 
momentum increase of 1ft per particle, (b) Energy levels of 
a 2D harmonic oscillator, (c) One of the possible spin-flip 
processes, bringing both particles to higher angular momen- 
tum states, (d) Eventually, after repeated application of the 
driving scheme, all particles occupy the lowest Landau level. 
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quency uj and an axial frequency ui z . For strong z con- 
finement hu z 3> E-p, where E-p is the Fermi energy de- 
rived below, the system is effectively 2D, see Fig. fa. 
The interactions between the particles are described by 
the dipolar interaction potential 



Vdd(r) 
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Tj is the relative distance between the two 



where r = r 

particles. Note that a weak s-wave scattering length does 
not change the general behavior of our transfer scheme 
and is ignored in the following. The dipole moment 
A* = (J-BgS/H = [i^gcrjl is given in terms of the Lande 
factor g and the Pauli matrices. By integrating out the 
fast motion perpendicular to the xy plane, taking the 
limit uj z — > oo, and using the spin raising and lowering 
operators <j ± = (<j x ± ia y )/2 the interaction reduces to 



V dd (r,cf>) 



a 
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where r, (j) are polar coordinates in the xy plane and 
Cdd = MoMbsVI^ 71 " characterizes the strength of the in- 
teraction. The dipolar interaction features three different 
processes. The first term proportional to cfcj describes 
spin-preserving collisions, while the second term cxj 
accounts for spin-exchange collisions. These terms con- 
serve separately the total spin and the total angular mo- 
mentum. Finally, the third operator e 2z ^ <J~a~ describes 
the relaxation process that transfers spin to orbital an- 
gular momentum, see Fig. fb,c. The sum L + S is still 
conserved and the spin flip leads to an orbital motion 
with an increase of relative angular momentum of 2h. 

It is this process that allows us to drive the dipolar 
particles to higher angular momentum states. Assuming 
the gas is initially in a spin-polarized state with the ex- 
ternal magnetic field pointing in the positive z direction, 
the particles will undergo spin relaxation when the field 
is adiabatically ramped through zero and finally pointing 
in the negative z direction. During this adiabatic ramp- 
ing, the total orbital angular momentum is increased by 
NH with iV the number of particles in the system. For 
the goal to reach the lowest Landau level regime, it is re- 
quired to transfer L* = N(N~ l)/2-h angular momentum 
to the orbital degrees of freedom, as described below. It 
is therefore necessary to reverse the magnetic field and 
the spins to their original position, in a way that guaran- 
tees repeated application of the transfer scheme without 
affecting the orbital angular momentum. 

To achieve this, we propose rotating the magnetic field 
by 180° around an arbitrary axis lying in the xy plane 
(say, the y axis), slowly enough such that the spins ro- 
tate adiabatically, but fast enough such that the orbital 
degrees of freedom cannot follow. To satisfy the adia- 
baticity with respect to the spins and diabaticity with 
respect to the external degrees of freedom, the speed of 
the rotation 7 rot has to satisfy lo <C 7 ro t "C WLj where 



= gunE/h is the Larmor frequency. After the ro- 
tation, the magnetic field has enclosed a D-shaped path 
in the xz plane. The spins are now pointing upward 
(in analogy to Fig. Ib but with increased angular mo- 
mentum) and the transfer scheme can be applied again. 
Multiple repetitions are realistic and only limited by the 
finite lifetime of the trapped ensemble. 

High angular momentum states are indeed related to 
the quantum Hall regime, as there is a close connection 
between the Landau levels and the states |n, m) of a 
two-dimensional harmonic oscillator in terms of a radial 
quantum number n = 0, 1, . . . and angular momentum 
hm, see Fig. fb. In particular, the ground state of N 
fermions filled into the harmonic oscillator with the con- 
straint L = L* is given by the many-body state 



N-l 



Here Zk = (x^ + iyk)/lno are complex coordinates of 
the particles, A is the antisymmetrizer, M is a nor- 
malization constant, and Zho = yh/mu is the har- 
monic oscillator length. This wave function is equiv- 
alent to the Laughlin wave function for integer filling 
v = 1, with Zho = \J h/muj replacing the magnetic length 
V2l m = yj2hc/eB for electronic systems. Quite gener- 
ally, the states with n — and m > correspond to the 
states in the lowest Landau level, see Fig. Id. To reach 
the LLL regime, we have to repeat the transfer scheme 
at least L*/Nh = ( N - l)/2 times. 

To quantify a single transfer process, our first aim is 
to calculate the total energy of N harmonically trapped 
fermions for a fixed total angular momentum L (polar- 
ized state, one spin component). For the noninteracting 
system, the energy can be obtained by simple summa- 
tions. We start with the ground state for L = 0, where 
all energy shells up to the Fermi energy are completely 
filled. The energy of the single particle states \n,m) is 
given by E nm = hoj(2n+ |m| + 1). To avoid cluttering of 
notation, we introduce dimensionless quantities indicated 
by a ~ sign. These quantities are measured in oscillatory 
units. That is, energy in units of Hu>, angular momentum 
in units of H, lengths in units of Zho an d time in units of 
u) . The degeneracy of each energy level is simply given 
by g{E) = E. With N = J2d(E) = E F {E F + l)/2 the 
Fermi energy is determined by 



E, 



\(Vm + i-i)^V2N 



The total energy for N particles is then given by 



E(N) 



E 

E=l 
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which shows the known scaling of a trapped 2D Fermi 
gas [19]. To derive the total energy E(N,L) for L ^ 0, 
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we define N m as the number of particles with angular 
momentum m. The energy in terms of N m is given by 

E = J2 £ E nrn = ^N m (N m + \m\). (2) 

m n— m 

The exact ground state energy can be found combina- 
torially for small particle numbers by varying the N m 
for fixed N and L. The result for N — 10 is shown 
in Fig. 2. For larger particle numbers this method is not 
feasible, but an analytic solution can be found for large 
particle numbers. Then, we can treat N m as a contin- 
uous function. To find the minimum of (2) at fixed N 
and L, we introduce two Lagrange multipliers /i, Q for 
the conditions N = J2 m N m and L — J2 m N m m, respec- 
tively. Taking the functional derivative with respect to 
N m yields N m = (/} — \m\ + Om)/2. The parameters 
can be determined by solving the constraints and sum- 
ming from 7Ti— — —p,/(l + fl) to m + = /t/(l — f2), where 
N m± — 0. One finds 



J(2N) 3 + (3Lf ' y/(2N) 3 + (3L) 2 ' 

By using these relations and omitting correction terms of 
order l/L and y/N, we obtain the total energy 

E(N,L) = y{2N) 3 + (3L) 2 . (3) 

This result agrees with the exact behavior for L = as 
derived in (1), and even for particle numbers as small as 
N = 10 it is close to the exact ground state energy, as 
shown in Fig. 2. For L > L*, the minimization problem 
becomes trivial as all particles occupy the lowest Landau 
level. The energy is exactly given by E = L, which is also 
obtained asymptotically from (3) in the limit L> N . 

It is now possible to quantify the link between our ap- 
proach (fixed angular momentum) and rapidly rotating 
systems (fixed rotation frequency) explicitly. Both are 
connected by a Legendre transform and we should in fact 
interpret the Lagrange multiplier £1 = ^ as the rotation 
frequency. In a harmonic trap, the system becomes un- 
stable if exceeds the value of the trap frequency u, as 
the harmonic confinement in the rotating frame is effec- 
tively given by uj — O. The angular momentum 

_ (2A) 3 / 2 Cl 
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has a singularity for Q = O/cj = 1 and large values of 
L can only be achieved by tuning close to the critical 
value. It is this precise control on the rotation frequency 
that so far prevented the experimental realization of the 
quantum Hall regime in harmonically trapped gases. In 
contrast, for the present situation, the system is always 
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FIG. 2. Exact ground state energy (dots) for N = 10 parti- 
cles at fixed angular momentum L, compared to the approxi- 
mate expression (solid line) as given in (3). For L > L* = 45/1, 
the energy increases linearly. Inset shows L as a function of 
the rotation frequency Q. in the analytic model. L diverges 
at the critical rotation frequency Q = uj, when the rotation 
exceeds the trap frequency. 

stable as ft(L) < uj for all L. An arbitrary orbital angu- 
lar momentum can be transferred to the system by the 
ramping scheme with high precision. 

Starting from expression (3) for the energy, we are 
now able to describe the transfer process in the adiabatic 
limit. Let be the number of particles in the spin-up 
state and iVj, = N — iV-j- the particles in the spin-down 
state. We describe both components separately and write 
the total energy as E(N^, L t ) + E(N±, L±)+A- N± where 
we have introduced the Zeeman energy shift A = fisgB 
(energy measured with respect to the energy of the lower 
Zeeman state). We assume that every particle even- 
tually takes part in the transfer process (adiabaticity) 
and consequently one quanta of angular momentum is 
transferred per particle. Starting from the nonrotat- 
ing state at L = 0, this imposes the transfer condition 
+ L_i = L = N^h. Adding this condition with an- 
other Lagrange multiplier, one can quantify the transfer 
process as a function of A, see Fig. 3a. Coming from 
high fields where A > Ep, the transfer starts right at 
the Fermi energy. Note that during the transfer, while 
E-p > A > — Ep, both components (t, 1) rotate in the 
same direction. Eventually all particles get transferred 
to the lower spin state and the total angular momentum 
equals L = Lj. = Nh. 

To justify the adiabaticity assumption above, we sim- 
ulate the transfer process for small systems of few parti- 
cles. We include all interactions mediated by Vdd(f, <f>), 
and assume, that the strength of the interaction Cdd = 
(Cdd/^Ho)/^- 1 <C 1 is weak compared to the Landau level 
splitting. Then, only a few excited states have to be taken 
into account. The system dynamics is described by 

H = J2 l E nm + A(i) S a J c\c :i + ^ ]T V ijkl c\c] Cl c k 

i ijkl 

where each of the indices ijkl of the fermionic oper- 
ators labels a set of quantum numbers (n, m, a) and 
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FIG. 3. (a) Description of the transfer in the analytical model with TV = 100 particles for decreasing energy splitting A between 
the two components f and J,. The transfer starts at A = Ef with particles continuously being transferred into the J. state as 
A is lowered to —Ef- Notice that during the transfer, both components rotate in the same direction. The crossing Nf = TVj_ 
is not precisely at A = due to the initial bias, (b) Full simulation of the transfer scheme for TV = 4, 6, and 8 particles in 
the adiabatic limit 7 — > 0. As the Zeeman splitting A is tuned through zero, the angular momentum increases in steps of 2h, 
indicating the transfer of two particle at a time. The interaction strength is given by Cdd = 0.1. (c) Angular momentum at 
the end of the transfer for N — 6 particles at different values of the Landau-Zener parameter A = CcLj/7. The data points for 
different rates collapse onto a single curve. The solid line is a probabilistic model, fitted to the data points. 



A(t)/hu) = — jt is the time-dependent Zeeman shift, con- 
trolled by the linearly decreasing magnetic field. The 
calculation of the dipolar interaction matrix elements 
Vijki ~ Cdd is presented in the supplemental material. 
The only relevant parameters in this model are the trans- 
fer rate 7 = 7/0; and the interaction strength Cdd- For 
the perfect adiabatic transfer, in the limit 7 — > 0, we 
can find the instantaneous ground state of H as A de- 
creases. The results are shown in Fig. 3b for TV — 4,6, 
and 8 particles. The total angular momentum L(A) in- 
creases gradually from L = to L = Nh in steps of 2H, 
indicating that two particles are transferred at a time. 

To obtain results for a finite transfer rate 7, we sim- 
ulate the full time-dependent many-body problem. The 
total angular momentum L(t — > 00) at the end of the 
transfer for TV — 6 particles is shown in Fig. 3c for dif- 
ferent values of Cdd and 7. Remarkably, the data points 
collapse onto a single line using A = C^/7. This param- 
eter arises in the Landau-Zener formula of a single level 
crossing, and the collapse indicates that each pair trans- 
fer is dominated by an individual avoided level crossing. 
A simple model accounting for this behavior (solid line) 
describes the final angular momentum observed in the 
full simulation (see supplemental material). 

The preparation of the integer quantum Hall state with 
an orbital angular momentum of L* = TV(TV — l)/2 is fi- 
nally achieved by a sequence of ramping cycles: Starting 
with an unpolarized sample with the fermions equally 
distributed between the two spin states, i.e., TV-;- = iVj. = 
TV/2, a first transfer increases the orbital angular mo- 
mentum by only L = N/2. Then, N/2 — 1 subsequent 
cycles will transfer exactly the required orbital angular 
momentum to reach the integer quantum Hall state. 

In an experimental realization with 161 Dy atoms, the 
number of cycles can be significantly reduced due to the 
total spin of F = 21/2 in the hyperfine ground state. Al- 



though calculations for 22 internal levels are too complex, 
we expect no qualitative modifications, except that 21H 
of angular momentum are transferred per particle and cy- 
cle [20] . Two important experimental requirements are a 
precise magnetic field control [21] as well as a determin- 
istic preparation scheme for a certain particle number, 
as demonstrated in [22]. For the magnetic field ramp 
we can estimate an optimal minimum value for the rate 
7 = 2E-p/i c = 2\/2N /ujt c by observing that the Zee- 
man splitting has to be tuned at least once from E-p to 
— E-p within the experimental accessible time i c , which 
is limited by the lifetime of the atoms in the trap. The 
Landau-Zener parameter is finally given by 

, _ ut c ( Zddi \ 2 
~~ 2^2N V4feoy 

where the length Zddi — mfiQ^j, 2 /4irh 2 parametrizes the 
strength of the interaction [14]. In a setup with N ~ 10 
fermionic 161 Dy atoms, a long lifetime of t c = 10s and a 
radial frequency of u) = 3kHz are needed to reach values 
of A on the order of 1. We comment, however, that the 
transfer scheme works already for smaller values of A. 

A particularly interesting property of the integer quan- 
tum Hall state, potentially useful to detect the success- 
ful generation, is the perfectly flat density n — l/vr^ 
within a circular region of radial size \/NIho- In addi- 
tion, it is possible to reach states with L > L* by con- 
tinuing the transfer scheme. In this regime, highly corre- 
lated ground states appear that are closely connected to 
fractional quantum Hall states, see [12] for a discussion 
in the context of rotating systems. Consequently, the 
presented method of dipolar relaxation allows for the ex- 
ploration of integer and fractional quantum Hall states, 
but avoids the experimentally challenging requirement of 
precise control of the rotation frequency by directly tun- 
ing the orbital angular momentum. 
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Matrix elements of the dipolar interaction 

In this section, we show how to simplify the matrix el- 
ements of the dipolar interaction Vijki — (ij\Vdd\kl) with 
each index ijkl representing a set of quantum numbers 

i = (rii, m,i, <7j) and 
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the dipolar interaction in 2D in terms of the polar coor- 
dinates of the relative vector r = r*i — r 2 between two 
particles. The spin-part is easily resolved and therefore 
we can concentrate on matrix elements of the form 



r- 



n 1 m 1 n 2 m 2 



with Am = 0, ±2. It is useful to change the basis to 
center-of-mass and relative coordinate states with 



Q = (ri +r 2 )/V2, q=(r 1 -r 2 ) 



Notice the symmetric definition with the additional fac- 
tor of 1/V2 compared to the usual definition of the rel- 
ative vector r = \/2q. Due to the quadratic character 
of the potential, Q and q are subject to the same har- 
monic potential. Thus, the product state \n1m1n2m2) = 
\n\m\) \n2Tn2) can be decomposed in terms of harmonic 
oscillator states \NM) \nm) of the Q,q coordinates via 



|n 1 min 2 m 2 ; 



E 



N,M,n,m 



where the r" imi,l2m2 
wnere me 1 NMnm 



are called Talmi-Moshinsky coeffi- 
cients [1, 2]. Since the center-of-mass is not affected by 
the interaction, the relevant matrix elements are 



(V2q) 3 



nn: ) = S m+ Am,m' j dg 




where the radial functions i?™ are given in terms of the 
generalized Laguerre polynomials L ™ as 



The matrix element V Am is thus 

yAm _ ,rp*\n'i m 'i n ' "~ 

V ~ 2-^i >• >NMn'( 
N \M ,n,n' ,m 
00 



mnimiTi2m2 
(m+Am) ' 1 NMnm 



R?, +Am (q)R™(q)q 
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where the remaining integral can be calculated analyti- 
cally for specific n, n', m and Am. 



Probabilistic model for Ljv 

Here we derive the total amount of angular momen- 
tum after the transfer by assuming that each 2-particle 
process is described by an independent Landau-Zcncr 
(avoided) crossing. That is, we neglect any interference 
effects. For each Landau-Zener process, we define the 
probability to transfer the n-th pair of particles by P n — 
1 — e -A / A " with A = C* 2 d /7 the Landau-Zener parame- 
ter and A„ an effective coupling strength, describing the 
n-th pair-transfer process. The total angular momentum 
for N particles after one cycle is then given by weighting 
each possible outcome (L = 0, L = 2, . . . , L = N) by the 
respective probability 

N/2-1 

L N = 2nP 1 ---P n {l- P n+1 ) + NP 1 ---P N/2 

n=l 

= 2Pj (1 + P 2 (1 + P 3 (1 + • • • (1 + Pn/2)))) 
For N = 6 particles this reduces to 



U = 2P 1 (1 + P 2 (1 + P 3 )) 
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The assumption of independent crossings can now be 
justified a-posteriori. By fitting Lq to the simulation 
data we find Ax = 0.0056, A 2 = 0.025, A 3 = 1.74 with 
Ai <C A 2 <C A 3 . While we suspect this approximation to 
break down for larger N, the model describes the trans- 
fer for small particle numbers in good agreement with 
the simulation. 
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